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Abstract. A new version of Monte Carlo generator ELRADGEN for simulation of 
real photon emission in elastic electron-proton scattering is presented. The exten- 
sions in the new version include opportunity to deal with polarized particles: lon- 
gitudinally polarized electron and arbitrary polarized proton. Simulation strategy, 
specifications of used kinematics, structure of the contributions to the observed 
cross section, cross-checks, and numerical results for BLAST experimental setup 
are presented and briefly discussed. 

1 Introduction 

Exclusive real photon production in lepton-nucleon scattering plays a rather important role 
in the investigation of the nucleon structure. The measurement of this process in different 
kinematical regions allows researchers to obtain the information about the generalized parton 
distributions [l|2j and the generalized polarizabilities |3|4j . In some cases this process appears 
as background to lepton-nucleon scattering in elastic [5] and inelastic [3] channels. 

Here we present a new version of the Monte Carlo generator ELRADGEN. The previous one 
[7] was developed for simulation of real hard photon emission from the lepton legs as background 
to unpolarizcd elastic electron-proton scattering. In the new version 2.0 we extend the generator 
to deal with initial polarized particles: longitudinally polarized lepton and arbitrary polarized 
proton. Both new and previous versions of this generator are based on the results of ref. [8] 
which were obtained using the Bardin-Shumeiko covariant approach for the extraction and 
cancellation of the infrared divergence [9] . 



2 Kinematics and Method of Generation 

For the simulation of exclusive radiative events in polarized electron-proton scattering 

e(h, b)+p(pi, V)^e(k 2 )+p'(p 2 ) +1 (k) (1) 

(k 2 = 0, k\ = k 2 = m 2 , v\ — v\ — M 2 ) we choose three kinematic variables: a transfer 
momentum squared t = —{ki — k 2 — k) 2 , the inelasticity v = (p 2 + k) 2 — M 2 , and the azimuthal 
angle cj>k between the planes (q, k) and (ki,k 2 ) depicted in Fig. [1] (a). Together with the 
kinematic variables characterizing Born contribution to elastic scattering 

Q 2 = -q 2 = -(ki-k 2 ) 2 , S = 2k lPl , <f>, (2) 
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Fig. 1. Three-vectors decomposition in Lab system: a) momenta leptons and photon; b) polarized 
vector for longitudinally polarized lepton £l and arbitrary polarized nucleon rj. 



they represent a full set variables for reconstruction of the four-vectors of all final particles in 
any frame. 

The four-vectors £l and r\ in ([1]) describe the longitudinally polarized electron and the 
arbitrary polarized proton, respectively. The proton polarization is characterized by two angles 
Orj and (firj as presented in Fig.[T](b). The explicit expressions for polarized vectors can be found 

Simulation of exclusive radiative events requires a separation of the total (or the radiatively 
corrected) cross section cr fcJ3 into the contribution of real hard photon emission <J r ad(vmin) 
and remaining part containing Born, soft photon, and additional virtual particle contributions 
&bsv (v m i n ). It can be performed by introducing a separation parameter, namely, the minimum 
inelasticity value v m i n that can be associated with missing mass square resolution of the detector 
[TO] when the final proton are not detected. The sum of these two positive parts 

Vobs = Vrad{Vmin) + &BSv{ v min) (3) 

does not depend on v m i n while <JBSv{v m in) and o ra d(v m in) do. The numerical details of these 
dependencies are presented in the next section. 

The explicit expressions for these two contributions are similar to those for cr ra d(v m i n ) and 
<^non-rad{vmin) from [7] with two exceptions: i) abbreviation BVS is used instead of non — rad 
for the remaining part of the cross section and ii) two structure functions representing the 
contributions of the polarized parts of the cross sections are additionally used. 

The strategy for simulation of an event can be outlined as follows. Two contributions 
&rad(v m in) and cfBSvivmin) to <j b s are calculated using a predetermined value of the v m i n . 
Then the channel of scattering (i.e., the process with or without real hard photon emission) is 
simulated according to the partial contributions of these two parts to the radiatively corrected 
cross section. If the channel with the real hard photon emission is chosen, three photonic vari- 
ables t, v and <pk are simulated according to their calculated distributions as they contribute 
to <Jrad(v m in) (see Fig |2] and ref. [JJ for details). Finally using these variables together with 
Born ones ([2]) which can be simulated according to the Born cross section or be externally 
predetermined, the four-vectors of all final particles in any frame are reconstructed. 

The cross section of the process ([1]) when the real hard photon emitted from the lepton legs 
is expressed through the nucleon form factors which depend only on one integration variable, 
namely, on t. Therefore to have an convenient opportunity to apply this generator to different 
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Fig. 2. Histograms (points) and corresponding probability densities (solid lines) for variables describing 
the exclusive real hard photon production in polarized electron proton scattering at BLAST kinematic 
conditions 11 {E beam = 850 MeV, Q 2 = 0.2 GeV 2 , 0„ = 48°), with <p = fa, PlPn = -1 and 
v mm = 10" 2 GeV 2 . 



fits or models of nucleon form factors, the integration over t has to be used numerical only, while 
the analytical integration over the other photonic variables v and 4>k is possible and allows to 
speed up the process of event generation. The analytical integration over these two variables 
was used in previous version of this generator for unpolarized scattering \j[ . However when we 
deal with the arbitrary polarized proton the analytical integration over v is a rather difficult 
due to non-trivial dependence of transverse component of the proton polarized vector on this 
variable. Therefore, in present version of the generator the analytical integration over v is used 
for unpolarized particle scattering only. 



3 Numerical test 



One cross-check which has to be done first is to investigate how well the simulated distributions 
of photonic variables t, v and 4*k reproduce those theoretically calculated. Fig. [2] provides such 
an illustration for the polarized electron proton scattering at BLAST kinematic conditions llj: 
the longitudinally polarized electrons with energy 850 MeV scatter off the proton that polarized 
under angle 9 V — 48° to the lepton beam. Plots for the photonic variable distributions contribut- 
ing to <J ra d{v m i n ) demonstrate excellent quality of the simulation procedure. Here Pl (Pn) is 
the degree of electron (proton) polarization, the solid lines corresponds to probability densities 
calculated numerically and the points show the event distributions simulated by ELRADGEN. 
Nucleon form factors are taken from [12]. In the left plot the theoretical ^-distribution is pre- 
sented when the integration over two other photonic variables has been performed. The sharp 
peak corresponds the situation when t = Q 2 . On the second plot one can see ^-distributions at 
the fixed t and integrated 4>k variables. Here the peak at t = 0.2 GeV 2 near small v corresponds 
the infrared divergence (that cut off by v rn i n ). The peaks on other two distributions correspond 
the situation when the real photon emitted along the momentum of scattered lepton. The right 
plot presents ^-distribution for fully differential contribution of real hard photon emission. 
As one can see, the momenta of the most of emitted photons are concentrated near scattering 
plane i.e. when cfik = 0. 

Two studies were illustrated in Tab. [T] : i) the investigation of the w m i„-dependence of 
&rad(v m i n ) , a bsv {v m in) and their sum; and ii) the comparison of the ratio of the radiatively 
corrected cross section a b s to the Born contribution <ro obtained by our generator and Fortran 
code M ASCARAD [8J with the same input parameters as a simplest comparison of these two 
codes. Specifically, Tab.Q]demonstrates that the observable cross section does not almost change 
with decreasing v m i n from 1 to 10~ 4 GeV 2 while its components cr ra d(v m in) and <JBSv( v min) 
change essentially: <7 r ad(v m in) increases and <7BSv{vmin) decreases. The comparison with MAS- 
CARAD results in a good agreement as well. 
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Table 1. The ?; m m-dependence of the radiative, BSV and observable contributions to electron-proton 
scattering with polarized target for different spin orientation in the Born units and results of comparison 
with MASCARAD [8j at BLAST kinematic conditions [11] (E beam = 850 MeV, Q 2 = 0.2 GeV 2 , 
n = 48°), with <t> = 4> n . 
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4 Conclusion and Outlook 

In the present report the new version of the Monte Carlo generator ELRADGEN for simulation 
of real photon events within elastic electron-proton scattering generalized for longitudinally 
polarized lepton and arbitrary polarized target is presented. 

Numerical test of new version of this code shows a good agreement with the Fortran code 
MASCARAD [8] and reveals lack of dependence on minimum inelasticity value v m i n with 
accuracy up to 1%. Besides we found that the distributions of generated radiative events are 
in coincidence with corresponding probability density. 

The present approach is rather general and can be extended in many other different ways 
including i) the development of this generator for transferred polarization from lepton beam to 
recoil proton [13j for measurement of electromagnetic form-factors of the proton in polarized 
scattering [14|15j ; ii) its further generalization for the investigation of electroweak effects such 
as axial form factors of the nucleon [16] and parity violation elastic scattering [17] ; and iii) its 
generalization for practical involvement in the experiments with the measurement of generalized 
parton distribution [l|2j as well as generalized polarizabilities [3|4j . 

We would like to acknowledge useful discussion with E. A. Kuraev. The one of us (A. I.) would like 
to thank the staff of MIT Bates Center as well as O. F. Filoti for their generous hospitality during his 
visit. 
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